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Abstract: We study classical production of soliton-antisoliton pairs from colliding 
wave packets in (1 + l)-dimensional scalar field model. Wave packets represent mul- 
tiparticle states in quantum theory; we characterize them by energy E and particle 
number A^. Sampling stochastically over the forms of wave packets, we find the en- 
tire region in {E, N) plane which corresponds to classical creation of soliton pairs. 
Particle number is parametrically large within this region meaning that the proba- 
bility of soliton-antisoliton pair production in few-particle collisions is exponentially 
suppressed. 
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1. Introduction 

Wonders related to classical dynamics of solitons in non-integrable models surprised 
theorists for decades [1-4]. Intriguing long-living bound states of solitons and anti- 
solitons — oscillons — are found in a variety of models [5-9]. Another interesting 
example is kink-antikink annihilation in (l + l)-dimensional 0"^ theory which displays 
chaotic behavior [10-12]. Most of these phenomena are explained qualitatively by 
reducing the infinite number of degrees of freedom in field theory to a few collective 
coordinates [13-15]. Then, mechanical motion along the collective coordinates shows 
whether soliton evolution is regular or chaotic. 

Recently [16, 17] a question of kink-antikink pair production in classical 
wavepacket scattering was addressed^, cf. Refs. [19-22]. The interest to this question 
stems from the fact that, within the semiclassical approach, wave packets describe 
multiparticle states in quantum theory. Studying kink-antikink creation from wave 
packets one learns a lot about the quantum counterpart process: production of non- 
perturbative kink states in multiparticle collisions. The prospect of Refs. [16,17] was 
to describe a class of multiparticle states leading to classical formation of kinks. 

^Related yet different problem is creation of kink-antikink pairs from wave packets in the back- 
ground of preexisting kink [18]. 
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Due to essential nonlinearity of classical field equations the process of kink- 
antikink creation cannot be described analytically^ and one has to rely on numerical 
methods. A difficulty, however, is related to the space of initial Cauchy data which is 
infinite-dimensional in field theory. Because of this difficulty the analysis of Refs. [16, 
17] was limited to a few-parametric families of initial data. 

In this paper we explore the entire space of classical solutions describing soliton- 
antisoliton pair creation from wave packets. To this end we sample stochastically 
over the sets of Cauchy data and obtain large ensemble of solutions, cf. Refs. [23,24]. 
We select solutions evolving between free wave packets and sohton-antisoliton pair 
and compute the energies E and particle numbers N of the respective initial states. 
In this way we obtain the region in {E, N) plane corresponding to classical creation 
of solitons. We are particularly interested in solutions from this region with the 
smallest A^. 

The model we consider is somewhat different from the standard 0"^ theory used 
in Refs. [16,17]. We do study evolution of a scalar field in (1 + 1) dimensions but 
choose nonstandard potential V{(f)) shown in Fig. |l]a, solid line. The reason for 




Figure 1: (a) Potential V {(/)). (b) Soliton and its spectrum. 

the unusual choice is chaos in kink-antikink scattering in c/)^ theory: it would be a 
venture to try applying new method in a potentially nontrivial chaotic model. We 
will comment on generalizations of our technique in the Discussion section. 

With the above set of classical solutions we test the method of Ref. [25] where 
classically forbidden production of kinklike solitons in the same model was studied. 
Namely, we compare the boundary of the "classically allowed" region in {E, N) plane 

^In particular, collective coordinates cannot be introduced since solitons are absent in the be- 
ginning of the process. 
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with the same boundary obtained in Ref. [25] from the classically forbidden side. 
Coincidence of the two results justifies both calculations. 

The paper is organized as follows. We introduce the model in Sec. § and explain 
the stochastic sampling technique in Sec. ^. In Sec. § we present numerical results 
which confirm, in particular, results of Ref. [25]. We conclude and generalize in 
Sec. |. 



2. The model 



The action of the model is 



S = ^ I dtdx [{d^,(l)f/2-V 
9 t/ 



(2.1) 



where 0(t, x) is the scalar field; semiclassical parameter g does not enter the classical 
field equation 

{d^ -dl)<P = -dV{<P)/d<P . (2.2) 

We assume that the potential V{(f)) has a pair of degenerate minima (/>_ and 0+. Then 
there exists a static solution of Eq. ( |2.2| ) — topological kinklike soliton 05 (x) shown 
in Fig. |I]b. Antisoliton solution 4>a{x) is obtained from (ps{x) by spatial reflection, 

(f)A{x) = (psi-x). 

We consider classical evolutions of 0(t, x) between free wave packets in the vac- 
uum 0_ and configurations containing soliton-antisoliton pair. Initial and final states 
of the process are shown schematically in Fig. |[ We restrict attention to P 
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Figure 2: Classical formation of soliton-antisoliton pair. 



symmetric solutions, (f){t,x) = (j){t, —x). This is natural since soliton and antisoliton 
are symmetric with respect to each other. 

In what follows we solve Eq. ( p.2|) numerically. To this end we introduce a 
uniform spatial lattice {xj}, i = —N^, . . . of extent —L^ < Xj < L^. At lattice 
edges X = ±Lx we impose energy-conserving Neumann boundary conditions (9^0 = 0. 
We also introduce a uniform time step At. Typically, = 15, A^^; = 400, At = 0.03. 
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Discretization of Eq. ( p.2| ) is standard second-order^. We take advantage of the 
reflection symmetry x — )■ —x and use only one half of the spatial lattice. 
We consider the potential 

r(0) = i(0 + l)2 

where dimensionless units are introduced; W{x) = e~^^{x + x^ + x^), a = 0.4. The 
value of V is chosen to equate the energy densities of the vacua, v ~ 0.75. The 
potential ( p.3|) is depicted in Fig. ^a, solid line. We denote the masses of linear 
excitations in the vacua 0_ and (^+ by m_ and m+, respectively. 

Our choice of the potential is motivated in two ways. First, we have already 
mentioned that kink dynamics in the standard (p^ theory is chaotic [10-12]. The 
source of chaos hides in the spectrum of linear perturbations around the 0^ kink. 
The latter contains two localized modes: zero mode due to spatial translations and 
first excited mode representing kink periodic pulsations. Localized modes accumulate 
energy during kink evolution which is thus described by two collective coordinates. 
Mechanical model for these coordinates is chaotic [12], just like the majority of two- 
dimensional mechanical models. 

We get rid of the chaos by choosing the potential ( p. 3D where the spectrum of 
linear perturbations around the soliton contains only one localized mode. Due to this 
property soliton motion is described by one-dimensional mechanical system which 
cannot be chaotic. 

Let us compute the spectrum of the soliton in the model (|2.3|) . Consider small 
perturbations — (psi^) = S(f){x) ■ e^""^* in the background of the soliton. Equa- 
tion ( p.2|) implies, 

[-dl + U{x)] = u^5(l){x) , (2.4) 

where U{x) = V"{(f)six)) and nonlinear terms in 50 are neglected. Discretization 
turns the differential operator in Eq. ( |2.4| ) into a symmetric {2Nx + 1) x {2Nx + 1) 
matrix; we compute the eigenvalues {oj'jf''} of this matrix by the standard method 
of singular value decomposition. Several lower eigenvalues are shown in the inset in 
Fig. |l]b. One sees no localized modes between zero mode and continuum u^^^ > m_. 

Another, unrelated to the soliton spectrum, mechanism of chaos was proposed 
recently in Ref. [26]. This mechanism works under condition m+ < m_ which is not 
met in our model. 

The second reason for the choice ( |2.3| ) is linearization of classical solutions at 
large negative times. Interaction terms should be negligible in the initial part of 
the classical evolution; otherwise initial wave packets cannot be associated with the 
perturbative Fock states. However, (l + l)-dimensional solutions linearize slowly due 

^One changes to + — 2(j)i)/Ax^, where — (f>{xi). The time derivative d^cj) 

is discretized in the same way. 
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to wave dispersion. Brute force linearization would require large lattice which is a 
challenge for the numerical method. Our model is specifically designed to overcome 
this difficulty. At a <^ 1 the potential ( p.3|) is quadratic everywhere except for the 
small region 0^1. Wave packets move freely in this potential if their tops are away 
from = 1, see Fig. After collision the wave packets add up coherently and hit 
the interaction region 0^1. Below we find that all classical solutions of interest 
behave in the described way. We use a = 0.4 which is small enough to provide, for 
the chosen lattice size, linearization at the level of 1%. 

It is worth noting that the problem with slow linearization is absent in multidi- 
mensional theories because amplitudes of spherical waves in D > 2 decay at power 
laws with distance. In general (1+1) -dimensional model linearization can be achieved 
artificially by switching off the interaction terms of the potential at |x| > Lint- This 
corresponds to a physical setup where interaction takes place in a sample of length 

2-Z^mi- 

We compute the energy E and particle number of the initial wave packets in 
the following way. Since the wave packets move freely in the vacuum 0_, 

/2 f°° dk 
- / cos(kx) We-''^''^ + ate*'^'^*] as t ^ -oo , (2.5) 

TT Jq y/2Uk 

where we took into account the reflection symmetry and introduced the amplitudes 
ak] ul = k"^ + rn?_. Given the representation (pT5|), one calculates E and A^ by the 
standard formulas, 

2 2 
E=— dkujk\ak\'' , dk\ak\^ . (2.6) 

9 Jo 9 Jo 

Expression for A^ can be thought of as a sum of mode occupation numbers Uk = IflfcP, 
where the latter are defined as ratios of mode energies Uk\ak\'^ and energy quanta"^ 
Uk- Note that the energy E is conserved; it can be calculated at arbitrary moment 
of classical evolution as 

E = ^jdx [{dt<pf + (9.0)2 + 2V{<P)] . (2.7) 



In the case of free evolution this expression coincides with the first of Eqs. 
Needless to say that Eqs. ( |2.6|) can be used only in the linear regime; this is the 
practical reason for continuing solutions back in time until Eq. ( p.5|) holds. Below we 
check the linearity of classical solutions by comparing their exact and linear energies, 
Eqs. ( p.7|) and (p.6|). We characterize classical solutions by points in {E,N) plane. 

Expressions ( p.5|) , ( p.6|) are naturally generalized to the lattice system.^ One 
solves numerically the eigenvalue problem ( p.4| ), where U{x) = m?_, and finds the 



^In our units 1%= 1. 

^Discretization of Eq. (2.7) is standard second-order. 



- 5 - 



spectrum {6(j)k{x), u)k} of linear excitations above the vacuum 0_. In this way one ob- 
tains lattice analogs of the standing waves cos(A;a;) and frequencies = a/A;^ + m'^. 
Arbitrary linear evolution in the vacuum has the form 

(Pit, x) = 0_ + 5^ 6Mx) [ake-''^'' + 4e*"'=*] , (2.8) 

k 

cf. Eq. ( |2.5| ), where we used the eigenmode basis with normalization 

^ Ax50fc(a;i)5(/)fc'(xi) = , Ax = Xj+i - x^. (2.9) 

i 

One extracts the amplitudes from the classical solution 0(t, x) at large negative t 
by decomposing (pit, x), dt(pit, x) in the basis of S(pkix) and comparing the coefficients 
of decomposition with Eq. ( |2.8D . Summing up the energies and occupation numbers 
of different modes, one obtains 




(2.10) 



where Eq. ( |2.9| ) is taken into account. 

3. The method 

3.1 Modification of the potential 

It is difficult to select solutions containing soliton-antisoliton pairs in the infinite 
future. On the one hand, numerical methods do not allow us to extend (pit, x) all 
the way to t — +oo. On the other hand, soliton and antisoliton attract; taken at 
rest, they accelerate towards each other and annihilate classically into a collection 
of waves. Thus, we never can be sure that (pit, x) contains solitons at t — t- +oo, 
even if lumps similar to soliton-antisoliton pairs are present at finite times. We solve 
this difficulty by changing the value of v in Eq. ( |2.3| ) and thus adding small negative 
energy density (— 5p) to the vacuum 0+, see Fig. |l|a, dashed line. This turns soliton- 
antisoliton pair into a bubble of true vacuum 0+ inside the false vacuum (p_ [27-29]. 
Large bubbles expand at Sp > since attraction between the solitons in this case 
is surmounted by the constant pressure 6p inside the bubble. Thus, at 5p > we 
simply look whether solution (pit, x) contains large bubbles at finite t. In the end of 
calculation, however, we have to consider the limit 6p — )■ 0. 

We remark that solutions containing soliton-antisoliton pairs at t — +oo can 
be identified by other methods. Our way, besides being particularly simple, has the 
following advantage: at 5p > there exists a critical bubble [27] — unstable static 
solution (pcbix) lying on top of the potential barrier between the true and false vacua. 
Given the critical bubble, one easily constructs classical evolutions between the vacua. 
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Indeed, in the critical bubble attraction between the soliton and antisoliton is equal 
to repulsion due to 5p. Being perturbed, it either starts expanding or collapses 
forming a collection of waves in the vacuum 0_. Thus, adding small perturbation to 
the critical bubble and solving classical equations of motion forward and backward 
in time, one obtains the classical solutions of interest. Critical bubble at 5p = 0.4 is 
depicted in Fig. |]a. 




-10 -5 5 

t 



(a) (b) 

Figure 3: (a) Critical bubble and (b) solution </>(t, x) describing its classical decay; 
6p = 0.4. Red (dark gray) and blue (light gray) colors in Fig. ^b mark regions with (j) > 
and (j) < 0, respectively. 

Let us obtain a particular solution describing creation of expanding bubble from 
wave packets at 6p > 0. We solve numerically Eq. ( |2.4| ) with U{x) = V"{(j)cb{x)) 
and find the spectrum of linear perturbations {S(j)^i^''\x) around the critical 

bubble. This spectrum is shown in the inset in Fig. it contains precisely one 
negative mode 6(j)neg{x), u^^g < due to changes in the bubble size. The latter mode 
describes decay of the critical bubble, 

(f){t,x) ^ (pcbix) + BnegS(j)neg{x)sh{\uJneg\(t -to)) , (3.1) 

where we fix to = 0, B^eg = 10~^ in what follows. Using the configuration ( |3.1| ) 
and its time derivative as Cauchy data at t = 0, one solves numerically Eq. (|2.2| ) 
forward and backward in time and obtains (j){t,x), see^ Fig. The latter solution 

^Only the central parts of solutions are shown in this and other figures. 
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interpolates between free wave packets above the vacuum 0_ and expanding bubble. 
We compute the values of {g'^E,g'^N) ^ (6.1,4.4) by Eqs. ( f^.lOD and mark the 
respective point "cb" in Fig. 




4.4 



^ 4 



3.6 



.cb 



N . (E) 

FJE)=0 



5a 



5b 



10 



(b) 



Figure 4: (a) Stochastic sampling at r = 50, = 10^. (b) Lower boundary of the 
"classically allowed" region. In both cases 6p = 0.4. 

In Fig. we check the linearity of evolution at t — t- — oo by comparing the linear 
and exact energies of (j)(t,x), Eqs. (|2.10| ) and (|2.7|) . As expected, the linear energy 
coincides with the exact one at large negative times and departs from it when wave 
packets collide. In what follows we estimate the precision of linearization as fractions 
of percent. 

3.2 Stochastic sampling technique 

In this Section we sample over classical solutions with bubbles of true vacuum in 
the final state, see Refs. [23,24]. It is hard to pick up the initial Cauchy data for 
such solutions: most of the initial wave packets scatter trivially and do not produce 
expanding bubbles. Instead, we consider the data {0(0, x), (9^0(0, x)} at t = 0. We 
decompose these data in the basis of perturbations around the critical bubble, 

0(0, X) = (j)cb{x) + Aneg S(f)neg{x) + ^ Ak 6 (f)''^^\x) , (3.2) 

k 

dt4>{0, X) = Bneg \uJneg\ S4>neg{x) + ^ ^fc ^^'j^^^ ^<Pk''\^) ' 

k 

where the negative mode is treated separately. Note that any functions 0(0, x), 
dt(f){0,x) can be written in the form ( p.2|) . For each set of Cauchy data {Aneg, Ak, 
Bneg, Bk} we solve numerically Eq. ( p.2| ) and obtain classical solution 0(t, x). Due to 
instability of the critical bubble there is a good chance to obtain transition between 
the vacua 0_ and 0+. 
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Next, we study the region in {E, N) plane corresponding to classical formation 
of expanding bubbles from colliding wave packets. We are particularly interested 
in the lower boundary N = Nmin{E) of this region. Let us organize the artificial 
ensemble of solutions describing transitions between the vacua. The probability of 
finding each solution in our ensemble is proportional to 

p oc e-^^-^'^ , (3.3) 

where E and are the energy and initial particle number of the solution; r and 
are fixed numbers. At large positive 'i? solutions with the smallest dominate in 
the ensemble and we obtain the boundary Nmin{E) with good precision. Value of r 
controls the region of energies to be covered. 

We use Metropolis Monte Carlo algorithm to construct the ensemble ( |3.3| ). In 
our approach solutions are characterized by the coefficients in Eq. ( p.2|) ; condition 
^neg = is used to fix the time translation invariance of Eq. p.2| ). The algorithm 
starts from the solution ( |3.1| ) describing decay of the critical bubble; it has Bneg = 
10-2, Ak = Bk = 0. Denote the energy and particle number of this solution by 
{Eq,Nq). We pick up a random coefficient from the set {A^, B^, Bneg} and change 
it by a small step. The latter step is a Gauss-distributed random number with zero 
average and small dispersion'' a. Substituting the modified set of coefficients into 
Eqs. ( |3.2| ), we find 0(0, x) and (9^0(0, x). Then, solving numerically the classical field 
equation, we obtain the entire solution (j){t,x). We compute the values of {E,N) by 
Eqs. ( |2.10|) . Solution is rejected if it does not interpolate between the vacua 0_ and 



0_i_; otherwise we accept it with the probability 

precept = min (1, e-^^-'''^^) , (3.4) 

where {AE, AN) are differences between the new values of {E, N) and the values 
{Eq, Nq) for the solution we started with. If the new solution is accepted, we write it 
down and use its parameters {A^, B^, Bneg}, {E,N) — )■ {Eq,Nq) for the next cycle 
of iterations. After many cycles we obtain the ensemble ( ^.3] ) of accepted solutions. 

A typical run of the Monte Carlo algorithm is shown in Fig. ^ where the accepted 
solutions are marked by dots in [E, N) plane. The algorithm starts in the vicinity 
of the critical bubble, then moves to smaller A^ and finally arrives to the boundary 
Nmin{E) where the majority of solutions is found. 



4. Numerical results 

We perform Monte Carlo runs at different values of r and until the entire curve 
A^ = A^min(-E) is covered with solutions. In total we obtained 2 ■ 10^ solutions, where 
the value of r was ranging between and 10^; = 10^, 10^, 5 ■ 10^. 

'^We used a — 10~^, 10~^; the final result was insensitive to this number. 
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The boundary = Nmin{E) is constructed by breaking the energy range into 
small intervals Ai? = 0.01 and choosing solution with minimal N inside each interval. 
This is the result we are looking for: Nmin{E) gives the minimum number of par- 
ticles needed for classically allowed production of bubbles. It is plotted in Fig. ^3, 
sohd line. As expected, Nmin{E) starts from {E,N) = {Ecb,Ncb) and decreases 
monotonously with energy. At high energies N^iniE) is approximately constant. 
Note that the particle number is parametrically large in the "classically allowed" 
region, A^mm ~ This means, in particular, that the probability of producing 

the bubble from few-particle initial states is exponentially suppressed. 

Given the boundary NminiE), we check results of Ref. [25] where classically 
forbidden transitions between A^-particle states and states containing the bubble 
were considered. The probability of these processes is exponentially suppressed in 
the semiclassical parameter, 

Vn{E) ~ e-^^^^y^' , (4.1) 



where Fn{E) is suppression exponent. One expects that this exponent vanishes in the 
"classical" region A^ > Nmin{,E). We extract the boundary of the set Fn{E) = from 
the results of Ref. [25] and plot this boundary in Fig. (dashed line). It coincides 
with Nmin{E) within 0.5% accuracy; the agreement justifies both calculations. 

Let us look at solutions with almost-minimal initial particle number, A^ ^ 
Nmin{E). Two such solutions are plotted^ in Fig. |^, their parameters {E,N) are 
shown by circles in Fig. Hb. At t — )■ — oo the solutions describe free wave packets 




Figure 5: Solutions 0(t, x) at 6p 
(10, 3.8). 



0.4. (a) (g^E, g^N) ^ (7, 3.9), (b) (g^E, g^N) 



^Small ripples covering tlie solutions represent fluctuations due to stochastic sampling. 
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moving in the vacuum 0_. After collision the wave packets emit waves and form the 
bubble. 

The most surprising part of the evolutions in Fig. ^ is emission of waves during the 
bubble formation. One assumes that the role of these waves is simply to carry away 
the energy excess which is not required for the creation of bubble. Indeed, solutions 
at different E look alike, cf. Figs. |]a and Pd; besides, Nmin{E) is independent of 
energy at high values of the latter. In Refs. [30-32] it was assumed that there exists 
certain limiting energy Ei which is best for bubble creation. Then, classical solutions 
with minimal particle number a.t E > Ei are sums of two parts: non-trivial soft part 
describing bubble production ed. E = Ei and trivial hard part — waves propagating 
adiabatically in the soft background. Hard waves carry away the energy excess 
E — El without changing the initial particle number; this is achieved at small wave 
amplitudes and high frequencies, see Eqs. ( p.lO| ). 

Numerical results do not permit us to judge whether the limiting energy exists. 
We can, however, confirm the conjectured structure of high-energy solutions. Con- 
sider the energies = uJk\ak\^ / g"^ of the modes at t — )■ — oo, where the amplitudes 
and frequencies ujk are defined in Eq. ( p^) . In Fig. ^ we plot these energies for the two 
solutions depicted in Fig. ^ Soft parts of the graphs are almost coincident. Long 




Figure 6: Mode energies at t — )• — oo for the solutions in Fig. 

tail of excited high-frequency modes is seen, however, in the graph representing the 
high-energy solution ^3. The tail carries substantial energy while the corresponding 
modes propagate adiabatically and do not participate in nonlinear dynamics. This 
is precisely the behavior conjectured in Refs. [30-32]. 

It is reasonable to assume that solutions at arbitrarily high energies have the 
same "hard+soft" structure. Then, Nmin{,E) stays constant as — )■ +oo and classical 
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formation of bubbles is not possible at any energies unless the initial particle number 
is larger than Nmin{E = +00). 

Finally, we consider the limit 5p — )■ 0. A typical solution at small 6p is depicted 
in Fig. It describes creation of soliton and antisoliton which move away from 



Figure 7: (a) Solution at 6p = 0.02, {g^E.g^N) ^ (10, 4). (b) Lower boundary N.min{E) 
at 6p 0. 

each other at a constant speed. The boundaries Nmin{E) at 5p = 0.04, 0.02 are 
plotted in Fig. |^ (dashed lines). They are almost indistinguishable; thus, the limit 
6p exists. Extrapolating Nmin{E) to 6p = with linear functions, we obtain the 
region in {E, N) plane for classically allowed production of soliton-antisoliton pairs 
(above the solid line in Fig. Qo). All initial wave packets leading to classical creation 
of soliton pairs have the values of {E, N) within this region. The region in Fig. |^ is 
qualitatively similar to the regions at 6p > 0; in particular, Nmin{E) is constant at 
high energies. 



5. Discussion 

In this paper we studied multiparticle states leading to classically allowed production 
of soliton-antisoliton pairs in (l+l)-dimensional scalar field model. We characterized 
these states with two parameters — energy E and particle number A^; we have 
found the corresponding "classically allowed" region in {E, N) plane. There were 
two main ingredients in our technique. First, we added constant pressure Sp pulling 
soliton and antisoliton apart. This modification led to appearance of the critical 
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bubble — unstable static solution lying on the boundary between the 
perturbative states and soliton-antisoliton pair. Second, we applied 
stochastic sampling over Cauchy data in the background of the criti- 
cal bubble. We thus obtained large ensemble of classical solutions de- 
scribing formation of soliton-antisoliton pairs. Calculating the values 
of {E, N) for each solution we found the required "classically allowed" F, 
region. 

Our method is naturally generalized to higher-dimensional mod- 
els. For example, consider formation of t'Hooft-Polyakov monopole- 
antimonopole pairs in four-dimensional gauge theories [33,34]. Con- 
stant force dragging monopole and antimonopole apart is pro- 
vided [35] by external magnetic field H, see the figure. At H 7^ ~9mH 
there exists a direct analog of the critical bubble [36] : unstable static 
solution where the attractive force Fatt between the monopole and antimonopole is 
compensated by the external forces ±gmH {g^ is a magnetic charge of the monopole). 
One performs Monte Carlo simulation in the background of this static solution 
and obtains many classical evolutions between free wave packets and monopole- 
antimonopole pairs. 

A particularly interesting application of our technique might be found in the 
study of kink-ant ikink production in (l + l)-dimensional 0^ theory. In this model the 
boundary Nmin{E) is lowered [16,17] due to chaos. We do not expect any difficulties 
related to nontrivial dynamics of solitons. However, classification of initial data for 
kink-antikink formation may require modification of our method. 
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